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Abstract 



Inclusive production of Xcj in e + e annihilation is an excellent probe of the 



Q\ ■ role of higher Fock states in the production of heavy quarkonia. Within the 



non-relativistic QCD approach, contributions from the short-distance produc- 
tion of colour-octet cc pairs are significantly larger than those from colour- 
singlet production. At the same time, Xcj production rates are significantly 



^ ■ smaller than expected in the colour evaporation approach. Measurements of 

Xcj production at CLEO and future B-factories will thus constitute a major 

test of theoretical approaches to the production of heavy quarkonia. 
PACS numbers: 13.65.+L 12.38.Bx, 13.85.Ni 
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Typeset using REVTgX 



Calculations of heavy-quark production in e + e~ collisions are now approaching a level 
of high accuracy. Fifteen years after the calculation of the next-to-leading-order (NLO) 
perturbative QCD corrections to the total open heavy-quark cross section we have now 
seen the completion of the NLO corrections to three-jet cross sections with massive quarks 
0. Despite continuous efforts over nearly 20 years, such a precision has not yet been 
reached in the calculation of cross sections of heavy bound states f3HT2"ll . Predictive power 



depends on the understanding of the long-distance bound-state formation process, which 
is therefore at the heart of current quarkonium physics. In this letter we investigate the 
inclusive production of P-wave (orbital angular momentum L — 1) quarkonium states in 
e + e~ collisions. This is the first study of these processes since essentially |13|] all previous 
work was concerned with the S'-wave J/ip or T mesons. 

For a sufficiently large quark mass m, a quarkonium bound state of a heavy quark Q 
and its antiquark Q is a non-relativistic system. The spectroscopy of both charmonia and 
bottomonia is well described in non-relativistic potential models, where the quarkonia are 
considered to be pure QQ states bound by an instantaneous potential. The question arises of 
whether the production of heavy quarkonia is also dominated by this leading Fock state. Is 
the QQ pair produced, already at short distances, in a configuration that corresponds to the 
asymptotic valence Fock component of the quarkonium? If it is, then all non-perturbative 
information in the theoretical prediction reduces to a single number, namely the coordinate- 
space wave function at the origin, which can be extracted, e.g. from a potential-model 
calculation. 

It is nowadays widely believed that higher Fock states are of decisive importance. Yet 
there exist differing approaches when it comes to relating observable cross sections to short- 
distance QQ production amplitudes. In the colour evaporation model (CEM) |3[], the quarko- 
nium production cross section is prescribed to be a (process-independent) fraction of the QQ 
production cross section, below the physical threshold for the production of a pair of heavy- 
light mesons. In the non-relativistic QCD (NRQCD) approach jjHf , the transition from QQ 



to quarkonium is represented in terms of a multipole expansion, which leads to scaling rules 



for transition probabilities in terms of v, the mean heavy-quark velocity in the meson. Any 
given cross section is a linear combination of several non-perturbative factors, multiplied 
by process-dependent hard factors. To leading order in Aqcd/a* (/•* ~ mv ) the former are 
well-defined process- independent NRQCD matrix elements (MEs). 

We have calculated the e + e~ — > Xcj + X cross sections within the NRQCD approach, 
using values of non-perturbative transition probabilities determined from measurements of 
other reactions. On the one hand, the results strongly violate the process-independence of 
cross-section ratios assumed in the CEM. Thus already an experimental upper limit on Xcj 
production in e + e~ annihilation can establish that non-perturbative effects play a decisive 
role in bound-state formation. On the other hand, the NRQCD results are dominated by 
colour-octet production channels. Observation of the predicted shapes and normalizations 
of the cross sections will therefore provide striking evidence of the importance of higher Fock 
states and the scaling of their contributions with v. 

In NRQCD, any cross section is given as a series expansion in both a s (m) and v, to 
leading order in Aq CD //i. Let us first consider the case of S'-wave quarkonia, in particular 
the J/i[) (J PC = 1 ). The leading contribution in v is given by processes where the cc pair 
is produced in the leading, colour-singlet Fock state [T5|], |cci( 3 Si)): 



e + e -> cci( 3 5i) + cc , (1) 
e+e-^ccifSO+gg. (2) 

The order v 2 correction to J ftp production still involves only the leading Fock state and 
can be considered as a relativistic correction to the wave function. Subleading Fock states 
first contribute at relative order w 4 , where the short- distance production of a colour-octet 
cc pair is followed by dipole transitions to the final quarkonium. The chromo-electric (El) 
and chromo-magnetic (Ml) dipole transition scale as v 2 (El) and v 4 (Ml), at least for 



a s {m c v) / n <C 1 |TT|. Dominant colour-octet J/-0 production processes at e + e colliders are 

e + e" -> cc 8 ( 1 S'o) +g , (3) 
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e+e" -> C c 8 ( 3 Pj)+g , (4) 
e + e" -> cc 8 ( 3 Si) + qq (q = u, d, s) . (5) 

Processes ([||]) are enhanced by l/a s (m c ) relative to (JT[) , (fj), and (|5|) is enhanced, at large 
energies, by a logarithm hx(s/m%) with respect to ([!]) (which, in turn, dominates (0) by a 
power s/m 2 at large s). Current estimates of the long-distance MEs (O^g ( 2S+1 Lj)) | ]14| , 
which parametrize the transition from a cci j g( 25+1 -^j) state to the J/ip, are listed, for exam- 
ple, in [11 1 . At the centre-of-mass energy studied by the CLEO experiment at the Cornell 



Electron Storage Ring, y/s = 10.6 GeV, the cross sections for the colour-singlet processes 
(H) and (fj) are 0.20 pb and 0.35 pb, respectively, about one third of the experimental cross 



section o = 1.65 ± 0.25 ± 0.33 pb [16]. Hence there is room for colour-octet contributions, 
which we estimate as 0.50 pb for (|3|), 0.60 pb for (|J), and 0.01 pb for (Q). 

Conclusions on the presence of colour-octet mechanisms would, however, be premature 
since estimates of the total J/ip cross section suffer from large uncertainties. First, the 
colour-singlet MEs are known to 50% at best, and a factor of 2 uncertainty is certainly 
not exaggerated for the colour-octet MEs. Secondly, there could be large perturbative 
corrections. And finally, truly relativistic corrections of order v 2 may also be large, cf. direct 
J/ijj production in hadronic collisions |I7| . 

Therefore it has been suggested |§ to study the energy distribution da/dz of the J/ijj, 
where z = 2Ej/^/y/s. The predicted large colour-octet processes (§) and (|) are concentrated 
near z = 1 and should thus dominate the upper end point of the z spectrum. A dramatic 
change should be visible, notably in the polar angular distribution of the J/ip w.r.t. the 
beam axis. Unfortunately, measurements at large z are plagued by a large background from 
radiative ip f production |I6|. 



Another potentially useful observable is the J/ip polarization, which is measurable via 



the angular distributions of its decay leptons. Partial calculations already exist p,[12|, but 
the measurement is definitely not easy. 

We propose the study of Xcj production as a means of investigating the importance of 
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colour-octet contributions in quarkonium production. Colour-octet mechanisms will leave 
distinct footprints in the total cross section, and in the energy and polar-angle distributions. 
Contrary to S-wave quarkonia, the J ++ Xcj mesons receive contributions from one of their 
higher Fock states |ccs( 3 5'i)g) already at leading order in v: the associated El transition 
suppresses the cross section by v 2 , but this suppression is compensated because an S'-wave 
operator scales as 1/v 2 relative to a P-wave operator. 
Colour- singlet production of Xcj proceeds through 

e + e" -> cci ( 3 Pj) + cc . (6) 

Note that the contribution from 

e + e"->cci( 3 Pj)+gg (7) 

is zero for pure photon exchange in the s channel. Also the Z-exchange contribution is 
negligible, at low energies because it is proportional to s/Mf, and at high energies because 
(0) is suppressed by m 2 /s relative to (J6|). 

Colour-octet contributions to the energy distribution of Xcj away from z = 1 arise from 
(§) and the processes 

e + e" -> cc 8 ( 3 S'i) + cc , (8) 
e+e-^ccsC^O+gg. (9) 

All three processes possess the same scaling in v and a s (m c ) as the colour- singlet process 
(H). Two further colour-octet processes dominate near the z — 1 end-point, namely (|3|) and 
(|). The process 

e+e-^cca^O+g (10) 

is negligible for reasons similar to (|^). 

The calculation of processes is standard but tedious. We calculate the distribu- 

tions in both the energy and the polar angle 9 of the Xcj, 
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(i " T S(z)[l + a(z)] . (11) 



dz d cos 9 



Details will be presented elsewhere. Here we comment on the relation of our results and 
previous calculations. 

At high energies, m 2 /s — > 0, the cross section for @ reduces to 

g = 2a(e + e- cc) -L (^( 3 P J )) D^ Pj) (z) , (12) 

where D^-lcc( 3 Pj) are the P ar t° n i c colour-singlet charm fragmentation functions for c — > 



cci ( Pj) + c. Our expressions for these functions agree with [18 , 19 



The fragmentation limit of (M) is more involved. Writing it as 



dcr _ 1 

— = 2cr(e e — > cc) — 
dz m 



2a(e + e" -> cc) — (Of\*S x )) A>(*) , (13) 



c 



we observe that Dg(z) agrees with (3.6) in [|I9| if we take fx 2 = z 2 (1 — z) s in the result of 
p| . Note that (H) is not the sum of two fragmentation processes, e + e — > cc followed by 
c - * Xj or c — > Xj an d e + e~ — > ccg followed by g — ► xj, as is the high-energy limit of (|J). 
We can decompose (p~3|) as 

S = ^ (°8 J ( 3 ^i)){2a cE [D^ c . c(3Si) (z) + R 8 (z) 



+ 



1 dz dcr cgg 



y=z/z 



dy 

where a C c, C cg is the cross section for e + e~ — > cc, ccg; the gluon fragmentation function is 
simply DfX C c(3s 1 )( z ) = ^oc s 5(l — z)/2A, and -Dc^cc( 3 5i) 1S obtained by appropriately removing 
the NRQCD ME and changing the colour factor in the colour-singlet c — > J/ip fragmentation 

,(8) 

c^cc( 3 5i) 



function given in ||. Also, -Dc-Le( 3 s ) coincides with the equal-flavour limit of the b 



bcg( 3 S'i) + c fragmentation function of [18|. However, from (|14]) it is clear that -Dc^lcc( 3 5i) 
is not the colour-octet part of the c — ► XcJ fragmentation function, as conjectured in ||18|| . 
Interference terms present in the ccg( 3 >S'i) + cc diagrams for the equal-flavour case (the 
non-zero remainder term Rg) forbid a simple rescaling of the colour-singlet result. 
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Our cross section for (|[) agrees with if we appropriately replace the colour-singlet 
NRQCD ME and overall colour factors by their colour-octet counterparts. Finally, our cross 
section for (|5|) reduces to a form that is analogous to the ^-integrated expression given in 

For our numerical results we use m c = 1.48 GeV, a s = 0.28, a em = 0.0075, and values 



for the NRQCD MEs as given in Two MEs have not yet been determined from any 



experiment. We estimate their magnitude by scaling other known MEs: 

(or j es ))*v^{or j ( 3 s 1 )) 
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2.8 x 10~ 4 GeV 3 , (15) 



3 (Oi^^St)) 

1.8 x 10~ 3 GeV 5 . (16) 



+ , n _ 3 



The total Xci,2 cross sections are dominated by colour-octet processes, as shown in Table |. 
The (infrared-finite) colour- singlet parts are two orders of magnitude smaller than the J/ip 
cross sections. Adding in the colour-octet contributions leads to a marked increase. Still, 
the rates are only about 1/20-1/10 of the J/ip production rate. This confirms that Xcj 
production in e + e~ annihilation is indeed very sensitive to the power counting of long- 
distance effects in quarkonium formation. The non-observation of Xc's at CLEO |16[ (for 
integrated luminosities of about 3.1 fb -1 on the T(4S) and about 1.6 fb -1 off the resonance) 



already causes problems |20] for the CEM, where the Xcj '■ J/4> ratio is predicted to be the 
same as was measured in fixed-target [I7| and collider |!T] hadropro duct ion, i.e. about 1 
(5/3) for J = 1 (J = 2). 

The rates in Table | are for direct production, i.e. excluding the production of Xcj i n the 
decay of other hadrons. B-meson decay contributions can be removed experimentally, e.g. 
by a cut on z at CLEO. There is also a contribution from ip{2S) decays, ij){2S) — >• 7 + Xcj- 
We estimate the Xci cross sections to be 11 fb for process (0), 19 fb for @, 8.3 fb for (§), 19 fb 
for (H]), and 0.64 fb for (|5]). In total we hence expect about 60 fb (10% less for Xcz), with 



about one half due to colour-octet mechanisms. We conclude that direct XcJ production 
dominates over indirect one. 

Figure [I] shows the energy distribution of direct XcJ production. A steep rise of the 
distribution at large z signals the importance of the e + e~ — > ccs ( 3 -Pj) +g process. Signatures 
for the other colour-octet mechanisms are also clearly visible, e.g. a Xci '■ Xc2 ratio close to 
1 at z < 1, as opposed to the ratio < 0.5 obtained from colour- singlet processes alone. 

With a separate measurement of x&,2 and with sufficient statistics to determine the 
double-differential cross section d 2 a/dzdcos6, the angular distribution parameter a(z) will 
also serve to identify NRQCD production channels. As shown in Fig. |2|, the colour-octet 
value of a(z) for Xci is significantly lower than the colour- singlet value, whereas for % c2 the 
colour-octet value (at large z) is significantly higher than the colour- singlet one. 

In summary, inclusive Xcj production is a powerful tool to establish the size of higher Fock 
state contributions in the formation of heavy bound states. Different theoretical approaches 
rely on different power-counting rules resulting in markedly different cross-section ratios. 
For example, the Xc2 '■ ratio is as low as 1/100 if only cc pairs in the leading Fock state 
contribute, while it is of order 1 if bound-state formation proceeds dominantly through 
(colour-singlet or -octet) S'-wave cc pairs, as it does in the CEM. The velocity-scaling rules 
of NRQCD yield a ratio of about 1/12, since the presence of colour-octet processes induced 
by higher Fock states is much more pronounced for Xcj production than for J/ip production. 
This is partly because colour-octet processes enter Xcj production at the same order as the 
colour-singlet ones, and partly because C-parity suppresses the process e + e~ — > ccgg, which 
dominates J /if) production. 

Signatures for colour-octet contributions do not only show up in a dramatic increase in 
total Xcj production rates. Differences are also clearly visible in the energy and angular 
distributions. The expected rates could already be visible at CLEO with current statistics, 
and will definitely be measured in the near future at CLEO and at B-factories. 

M.V. wishes to acknowledge financial support from Suomalainen Tiedeakatemia, 
Vaisalan rahasto. 
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TABLES 

TABLE I. Integrated cross sections and angular coefficient a for e + e~ — > Xci,2 + X at 
\/i = 10.58 GeV. The colour-octet Xc2 cross sections are a factor 5/3 larger than the corresponding 
Xci ones. The cut z > 0.693 serves to exclude Xcj's originating from B-meson decays. 

All z z> 0.693 





a[fb] 


a 


a[fb] 


a 


CC + CCl( 3 Pl) -> Xcl 


18.1 


0.44 


15.3 


0.49 


CC + CCl( 3 P 2 ) -> Xc2 


8.4 


0.10 


7.5 


0.09 


CC + CC 8 ( 3 S'i) -> Xcl 


6.1 


0.24 


4.0 


0.35 


qq + cc 8 ( 3 5i) -> Xci 


15.6 


0.26 


11.7 


0.34 


gg + CCs^S 1 !) -> Xcl 


5.5 


-0.03 


4.2 


-0.05 


g + cc 8 ( 1 S'o) -> Xci 


4.1 


1.00 


4.1 


1.00 


g + CC 8 ( 3 Pj) -f Xcl 


45.6 


0.64 


45.6 


0.64 


Total Xci 


95.0 


0.47 


84.8 


0.53 


Total Xc2 


136.5 


0.46 


123.4 


0.51 
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FIGURES 



. dG(Y f )/dz [fb] 




FIG. 1. Energy distribution of <r(e + e - — ► Xci,2 + X) at yfs = 10.58 GeV: total, colour-singlet, 
and colour-octet Xci contributions, and total and colour-singlet Xc2 contributions. The peak from 
the e + e~ — > ccg( l So, 3 -Pj) + g processes at large z has been smeared over a range of width ~ v 2 . 
In a physical process such a smearing follows from the energy transfer in the non-perturbative 
transition. 
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